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Abstract. In the almost Friedmann-Lemaitre model of 
the Universe, the density parameter, fi m , and the cosmo- 
logical constant, f^A, measure curvature. Several linearly 
degenerate relations between these two parameters have 
recently been measured. Here, large scale structure cor- 
relations at ~ 100 — 150 h^ 1 Mpc are found in the co- 
moving three-dimensional separations of redshift z ~ 2 
quasars. These function as a comoving standard rod of 
length i^gg w 130 ± 10 ft. -1 Mpc. A local maximum in 
the correlation function at ps Lj^gg/2 also appears to be 
significant. By combining separate radial and tangential 
standard ruler analyses, the lifting of the Q m — £1^ linear 
degeneracy within a single data set is demonstrated for 
the first time. 



Key words: cosmology: observations — cosmology: the- 
ory — distance scale — quasars: general — large-scale 
structure of Universe — reference systems 

1. Introduction 



In standard cosmology (Weinberg 1972), space is a 3- 
manifold flSchwarzschild 190C ; Luminet & Roukema 1999) 
of nearly constant curvature, i.e. space is approximately 
locally homogeneous. Geometric ways of measuring aver- 
age curvature include the use of phenomena of intrinsically 
fixed brightness or leng t h scale, i.e. of st andard candles 
( Pcrlmuttcr ct al. 1999 ;_ Ricss ct al. 199§|) and standar d 



rulers [e.g. |Mo et al. (1992; ); |Broadhurst fc Jaffc (1999|) 



Roukema Sz Mamon (2000 ); note also t he microwave back- 
ground angular s tatistical estimates ( Lange et al. 2000 ; 



Balbi et al. 2000 ) which can loosely speaking be thought 
of as 'theoretical' standard rulers], but have previously 
been found to lead to degeneracy in the fl m — plane 
(e.g. [Linewcaver 199S ) . However, inhomogeneities (pertur- 
bations) in density exist and can be statistically repre- 
sented by a Fourier power spectrum, and are believed to 
gravitationally collapse and form objects such as galax- 
ies and clusters of galaxies. Use of a characteristic feature 
of this spectrum at a scale 3> 10 hT 1 Mpc, the size of 
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Fig. 1. Projected comoving spatial distributions of the 
two quasar subsamples at redshifts 1.8 < z < 2.4, for 
(f2 m = 0.4, = 0.6). a. The right ascension subsam- 
ple, for a increasing from left to right, b. The declina- 
tion subsample for 5 increa sing from left to r ight. The 
discrete redshifts published ( lovino et al. 1996 ) are con- 
verted to continuous values by uniform random offsets. 
The latter are used below. Differences in magnitude lim- 
its, hence different number densities, are visible in the 
a subsample and are corrected for in both subsamples 
(cf Table @). What appears to be a cluster or a su- 
percluster in b. at (+200, 3150) ft. -1 Mpc is located at 
(a « l'W", S re -38.5°, z ps 2.24). 



the largest bound structures, should provide a comoving 
standard ruler for constraining the local geometrical pa- 
rameters (firru^A)- 
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Many observations of both galaxies and superclus- 
ters of galaxies indicate that the maximum in the power 
spectrum is peaked at a wavenumber where 
L L S S ± A^LSS re 130 ± 10 ft -1 Mpc ([Broadhurst ct aT 



Physics which could potentially be investigated in or- 
der to explain the feature at ^lss includes acoustic oscil- 
lations in the baryon-photon fluid before last scattering, 



1990j; paztanaga fe Baugh 199S|; [Einasto ct al. 1997b 



in high baryon density models (Eiscnstcin 1998; Meiksin 



Deng et al. 1996] ) (comoving length units). Since this stan- 
dard ruler should be valid independently of orientation, 
the different O m — degeneracies implied in the radial 
and tangential applications of the ruler should enable lift- 
ing of the degeneracy within a single data set, providing 
a potentially more powerful ruler than previous standard 
rulers or standard candles. 

It should be noted that while the existence of a broad 
maximum in the power spectrum is uncontroversial, not 
all observational analyses agree on whether or not there 
is a sharp feature in the power spectrum in this region in 
additi on to the broad maximum, and there is not yet any 



ct al. 1998; Peebles 1999), and features from Planck epoch 
physics which transfer to oscillations in the post-inflation 
power spectrum, for weakly coupled scalar field driven 
inflationary models ( Martin fc Brandenberger 2000[ ). At 
high redshift, the L L ss = 130 ±10 ft -1 Mpc feature [ "dis- 



tance" means comoving proper distance ( Weinberg 1972 ) 
t hroughout this paper] has been detected am ong quasars 
( Roukema fc Mamon 2000|; Deng ct al. 1994 ) and Lyman- 
break galaxies ( Broadhurst fc Jaffc 1999 ). 

Most applications of standard candles or standard 
rulers exploit either the radial redshift-distance relation 



maxim lum. For example, on one hand, lEinasto ct alj's et al. 200C|), but not both sim ultaneously. 



( Broadhurst fc Jaffe 1999) or the tangential redshift 
distance-angle relation (Pcrlmuttcr et al. 1999 ; Ricss et al 



clear a greement on the characteristic scale of the broad 1998[; |Roukema fc Mamon 2000| ; |Lange et al. 2000| ; |Balbi 



(1997) analysis of superclusters suggests a sharp peak at 
k = 2-k/L ~ 0.05 ft. Mpc -1 . But, on the other hand, while 
in the low redshi ft IRAS PSCz (point sou rce catalogue 
redshift) survey ( Sutherland et al. 1999 ), there is, at 
least, a broad maximum at around 0.02 ft Mpc -1 < k = 
2ir/L < 0.04 ft Mpc" 1 (fig. 1 of Sutherland et al. 1999) ), 
i.e. 320 ft -1 Mpc Z L > 160 ft -1 Mpc, there is no obvious 
sharp feature in this region. Nevertheless, there is a 2a 
significant sharp peak (see fig. 1 and comment in section 
6 of [Sutherland et al. 1999| ) which lies at 0.07 ft Mpc -1 < 
k < O.lOftMpc" 1 , i.e. 90ft 1 Mpc > L > 60 ft -1 Mpc in 
the PSCz. 

Possible reasons why the ^lss = 130 ft -1 Mpc feature 
found by other authors, if real, might have been missed in 
the Sutherland et al. (1999| ) analysis include 



Alcock & _Paczynski (19Y9) suggested the idea of using 
both constraints simultaneously, and suggested applying 
it at quasi-linear or non-linear scales, i.e. r < 10ft _1 Mpc, 
but did not discuss how to lift the degeneracy in the 
two curvature parameters (f2 m ,f2A) which remains af- 
ter using the local isotropy constraint, though they did 
suggest a theoretical method for separating out some of 
the peculiar velocity effects which are important at these 



small scales. 


Phillipps (1994 


)• 


Matsubara & Suto (1996), 


Ballinger, Peacock & Heavens (1996 


) and 


Popowski et al. 


(1998) followed up this idea, demonstrating 


; specific formu- 



(i) use of a different population (lacking in early type 
galaxies) , 

(ii) use of too large a bin size, i.e. too much smoothing, 

(iii) redshift distortion (velocity dispersion at small scales, 
smooth infall at larger scales) , though this is d iscussed 



briefly in the para graph preceding Sect. 3.1 of Suthcr 



land ct al. (1999| ) 



(iv) assumption of zero cosmological constant [e.g. if 
(f2 m = 0.3, = 0.7), then the length scale to 
cz = 45000 km/s is underestimated by 7%], 

(v) assumption of zero curvature by use of a power spec- 
trum analysis. 



lae and calculations regarding quasar pairs and the two- 
point auto-correlation functions of galaxies and quasars, 
including separation of local isotropy ('sphericity') and 
some of the peculiar velocity effects. 

However, by using a standard ruler in the linear 
regime, i.e. by using a feature at ^lss ~ 130 ft -1 Mpc, 
peculiar velocity effects become negligible, and the in- 
ability of this scale to evolve in a Hubble time provides 
an additional constraint in the (f^m, ^a) plane. For the 
r < 10 ft -1 Mpc auto-correlation function, the peculiar 
velocity effects are certainly important, and evolution in 
the length scale must be contended with, for example by 
model-dependent assumptions. 

2. Observational Analysis and Discussion 



In a previous analysis ( Roukema fc Mamon 200(i|) of a 



deep, dense, homogeneous quasar survey (Iovino et al 



The method of Roukema & Mamon (2000) did not as- 
sume the ^LSS feature to be sharp, though a broad feature 
would obviously have given a less significant (or maybe 
an insignificant) signal. In the present analysis, a reason- 
ably sharp (but low amplitude) feature consistent with 
^LSS ~ 130 ft. -1 Mpc is found. A secondary f eature con- 
sistent with the k ~ 0.08 ft -1 Mpc feature of Sutherland 



et al. (1999 ) is also found, but is not studied in detail. 



1996), only the tangential relation was used, to ensure that 
observational selection effects well known to cause non- 
cosmological periodicities in redshifts ( [Scott 199l] ) could 
not bias the result. In the present analysis of the high 
grade quasar candidate catalogue (Table El Fig. [I]), the 
technique of "redshift scrambling" (see Fig. ^| caption) is 
used to enable use of three-dimensional information in a 
way that avoids redshift selection effects. Since the red- 
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Table 1. The two angular subsamples of the quasar can- 
didate catalogue (Iovino ct al. 1996), defined by B1950 
limits (cvi < a < 0,2, 8\ < 5 < 82) and redshift limits 
1.8 < z < 2.4. In §2.1, these are further divided into red- 
shift subsamples 1.8 < z < 2.1 and 2.1 < z < 2.4, yield- 
ing the four independent subsamples whose analyses are 
shown in Figs |, | and |. In ^ (Fig. |), 1.8 < z < 2.4 is 
used, i.e. yielding just two subsamples. When generating 
Poisson simulations for a given subsample, further subdi- 
visions, at a' for the right ascension subsamples, and at 
8' ,8" for the declination subsamples, are defined in order 
to allow for variation in magnitude zero points or mag- 
nitude cutoffs between the different plates. That is, for a 
given subsample with several subdivisions i defined by a' 
or 5', 5", there are nt quasars in each subdivision i. For 
each subdivision i, m uniform random a or 8 values are 
generated in the appropriate a or 8 range for that subdivi- 
sion. This mimicks the magnitude limits. The subdivisions 
a' , 8' and 8" are then ignored when calculating £. Num- 
bers of objects N are indicated. 
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Fig. 2. Spatial two-point auto-correlation function £(r), 
for separations r in comoving units and (f2 m = 0.4, = 
0.6). The four angular /redshift subsamples are shown as 
dashed (5,1.8 < z < 2.1), dashed-dotted (6,2.1 < z < 
2.4), dotted (a, 1.8 < z < 2.1) and dashed-triple-dotted 
(a, 2.1 < z < 2.4) curves. The mean (£) and the stan- 
dard error in the mean o~/£\ are shown by the thick 
and thin solid lines respectively. The correlation func- 
tions are calculated in three-dimensional curved space via 
f(r) = {DD - 2DR/n + RR/n 2 )/ (RR/n 2 ) where DD, 
DR and RR indicate numbers of data-data, data-random 
and random-random quasar pairs respectively (Landy & 



shifts JiscH in the random and observational catalogues 
consist of exactly the same set of numbers, any redshift 
selection effects, which are independent of angle, should 
statistically cancel out ( psmer 1981 ) in calculation of the 
correlation function £(r) ( proth fc Peebles 1977 ). Since 
some real signal could also cancel out, in principle, this 
implies a conservative estimate of £(r), i.e. a lower limit 
to |C(r)|. 

2.1. Local maxima in the 3- dimensional correlation 
function 

Figs H and || show that, for reasonable values of (f2 m , f2 A ), 
a local maximum in the correlation function consistent 
with ij^gg = 130 ± 10 hr 1 Mpc is clearly present. By con- 



Szalay 1993), and n — 20 times more random points than 
data points are used. The random catalogues use (i) uni- 
form probability distributions in the two angular direc- 
tions (Table |^), and (ii) random permutations ("z scram- 
bles", Illb in |0 smcr 198l| ) of the observational set of 
redshifts. to avoid biases from redshift selection effects 



(Scott 1991). Bin size is 5/z." 1 Mpc and £ is smoothed 
by a Gaussian with a = 10 h^ 1 Mpc. The low values of £ 
at r < 20/i _1 Mpc are related to redshift roundoff error. 



( Groth fc Peebles 1977 ), which most often increases the 
precorrected values of £. Making an integral constraint 
correction usually requires assumptions on the intrinsic 
shape of £. To avoid these assumptions, it is more prudent 



trast, an (fi m — 1,^2 A — 0) 
local maximum Lu ulciii ar 



uin^eise wu 
L L aj 100 /t' 1 M 



-tiris- 



icquiiu 
111 cunlia- 



just to quantify the peak as a local maximum (Deng et al 



diction with the low redshift estimates of L^g g. A cor 
relation function consistent with the standard (Groth & 
Peebles 1977) galaxy-galaxy correlation function £(r) 
(r/5 h^ 1 Mpc) -1 - 8 is also present for r < 40/i _1 Mpc. 

What is the significance of the ^Lgg peak? This de- 
pends on where the zero level of correlation lies. In correla- 
tion function estimates where both sample and correlation 
are small, the problem of only having a finite volume of- 
ten requires a correction known as the integral constraint 



1994; Roukema & Mamon 200C). For a maximum at r ma x 
consistent with a peak at /L^gg, where / = 1, define £ r , 
£_ and £ + as the maximum value and the first minima be- 
low and above r max respectively, and take the maximum 
value of <j<£) (r) for r € [r max - / AL LSS , r max + / AL LSS ] , 
where AL L gg = 10/i _1 Mpc. Then, 



(S/N)j 



r-(£-+£+)/2 



max{cr (?) (r)} 
yields (S/iV)i = 3.2 for (f2 m ,fi A ) = (0.4,0.6). 
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Fig. 3. Correlation function as for Fig. g, but for a 
hyperbolic universe (f2 m = 0.3, 51a = 0.0) above, and a 
zero cosmological constant, flat universe (51 m = 1.0, fl\ = 
0.0) below. 




Fig. 4. Confidence intervals for rejecting the presence 
of a peak at -^Lgg ± AI/Lgg f° r various hypotheses on 
(fi m ,f2A). Rejection levels are 1 — P >0% (white), 1 — 
P >68% (hatched), l-P >95% (light cross-hatched) and 
f - P >99.7% (heavy cross-hatched). The l-P = 68% 
contour for testing a peak at L;Lgg/2 i test 
is shown in bold. For each pair (f2 m ,f2A), the peak po- 
sition is estimated as the value r max for which £(r) is 
maximum in 100 < r < 300 hr 1 Mpc. (For the ^LSS/ 2 
peak, 40 < r < 100 hr 1 Mpc is searched). The mea- 
surement uncertainty in the /£Lgg P ea k is estimated 
as Ar max = fAL LSS (for / = 0.5,1). The probabil- 
ity of finding r max close to /LLgg assumes Gaussian er- 
rors, i.e. P f ee erfc[|r max - fL LSS \/(V2a)], where a 2 = 
(Ar max ) 2 + (/AL LSS ) 2 . 



constraints for clusters and groups of galaxies, though to 
slightly higher Q m values than were fou nd in the purely 



Table 2. Probability (eq. ||) of finding one, two or three 
peaks among 200 random simulations at /-^Lgg where / = 
0.5, 1 and/or 2, for (n m , A ) = (0.4, 0.6). 



tangential analysis of the present survey ( Roukema fc Ma- 
mon 200C| ). the coincidence would be surprising if it were 



D sim D sim D sim 

J n R r\ l; 



psim psm 

0.025 0.01 0.335 < 0.005 0.005 0.005 < 0.005 



D sim 

JO. 5, 2 



D sim 

■M).5,l,2 



Fig. ^ shows that an automatic search for this peak, 
using a simple and robust method, i.e. using the value of 
r for which £ (r) is maximum over a very large interval in 
r, yields an approximately linear confidence band in the 
(Qui, £Ia ) plane. Since this band is consist ent with kine- 
matical ( |Carlberg et al. 1997 ; Mamon 1993) and baryonic 
fraction flWhitc ct al. 1993| ; [Hcnriksc-n fc Mamon 1994[ ) 



due to noise or systematic effects. 

Moreover, what appear to be peaks at L^gg/2 ± 
AL^gg/2 and near 2LLgg ± 2ALLgg are present, though 
to lower significance, with (S/N)o.5 = 2.3 and (S/N)2 = 
1.2 respectively, for (f2m,^A) = (0.4,0.6). Could any of 
the three peaks be induced by noise which has common 
statistical properties among all the four subsets, either due 
to shot noise or selection effects? Redshift selection effects 
have been removed by the use of z-scrambling. Angular 
selection effects may be present at a small level (see sec- 
tion 3.4 of Roukema & Mamon 2000), but are more likely 
to decrease the amplitude of any signal rather than in- 
troduce false correlations which mimic the signal found. 
Moreover, the converge nce of the separate tangential and 
radial analyses below (§ 2.2) suggest that the effects of an- 
gular selection are weak. 
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To test the properties of shot noise, random simula- 
tions were performed as before but substituted for the 
data. The probabilities that maxima can occur as close 
to and of at least the same signal-to-noise ratio as the 
observed values can be defined 



P 



f 



= P 



/ L LSSl < l r obs-/ L LSsl and 



(S/N)f m > (S/N)y 



obs 



(2) 



where 'sim' and 'obs' indicate simulations and observa- 
tions respectively, and the intervals in r are as above. 
Since one might suspect a single noise feature to cause, 
say, simultaneous features at two or three of the posi- 
tions, mimicking the observed signal, the probabilities of 
finding the peaks might not be independent, e.g. ~ 
P|im ~ psim > psimpsim for f ± j 2 g {0.5,1,2} is in 
principle possible. The results (Table ||) show that the 
hypotheses of either of the LLgg/2 or the -&LSS P e& ks oc- 
curring by chance are each rejected to 1 — p sim > 97%, 
and that of both occurring simultaneously is rejected to 
1 - P sim > 99.5%. The hypothesis of the the 2L LSS peak 
occurring by shot noise cannot be significantly rejected. 
These values vary throughout the (fi m , plane. 

Fig. | shows that the 68% > 1 - P > 0% confidence 
intervals for the L^gg/2 and £*Lgg peaks are consistent. 

2.2. Tangential versus radial correlations 

A standard ruler should not depend on orientation. Can 
use of both radial and tangential information lift the de- 
generacy of (f2m, ^a) constraints? To illustrate this, the 
full redshift interval 1.8 < z < 2.4 is used, but only pairs 
of objects within 30° of either the radial or tangential di- 
rections respectively are used. For the geometry of this 
survey, about 10% of pairs are radial and 60% of pairs are 
tangential according to this criterion. 

Fig. H clearly shows, both observationally and theoret- 
ically, the difference in the slopes of the radial and tangen- 
tial constraints at z « 2. A hyperbolic (f2 m + Qa — 1 < 0) 
universe is suggested by the 68% confidence limit, though 
a flat universe with fl\ = 1 — O m = 0.5 is only rejected to 
1 — P w 80% confidence, i.e. not significantly. The partially 
lifted degeneracy can be represented (at 68% confidence) 
as 

n m = (0.30 ±0.04)O A + (0.30 ±0.11), ft A <0.5. (3) 

This suggests a somewhat higher matter density and lower 
cosmological constant than other recent results. 

3. Conclusion 

The confirmation of the ^Lgg peak and the partial lifting 
of the ri m — degeneracy show that ongoing and fu- 
ture large quasa r surveys [in partic ular the 2 Degree Field 
Quasar Survey ( Boyle et al. 2000| ) and the Sloan Digital 




Fig. 5. Partial lifting of the £l m — Ha degeneracy. Confi- 
dence intervals are as for Fig. | for a. the radial (l-Pf 8 * 1 ), 
b. tangential (1 - P^ n ) and c. combined (1 - pfadptan) 
constraints, d. A simple model, for which the 95% > 
1 — P > 68% contours for radial (straight contours) and 
tangential (curved, nearly vertical contours) constraints 
showing near linear degeneracy are shaded as before 
(hatched regions); and the combined 68% > 1 — P > 0% 
contour is the half-ellipse-like shape nested between these 
(heavy contour), a.-c. The radial and tangential tests are 
performed for the ^Lgg ± AL^gg peak as in Fig. |[ ex- 
cept that only pairs oriented within 30° of the radial and 
tangential directions (respectively) are included in calcu- 
lation of £(r). Both angular subsamples over 1.8 < z < 2.4 
are used. The radial and tangential tests are assumed to 
be statistically independent, d. For the model, the red- 
shift interval Az and the angle 9 corresponding to ^Lgg, 
assuming that (fi m = 0.3, Ha = 0.0) and that z — 2.1, are 
calculated. For each pair (f2 m ,f2A)) the radial and tan- 
gential distance intervals implied by Az and 9 are calcu- 
lated, ignoring the initial assumption about Q m and 
These are treated as two independent 'experiments', and 
Gaussian probabilities p raQl 5 ptan f observing these val- 
ues, given £*LSS' ^-^LSS' ^ m and fl\ are calculated as 
before. The combined rejection is 1 — p rad ptan 



Sky Survey quasar sample (e.g. Fan et al. 2000| )] will have a 
much more powerful tool for local geometrical constraints 
than was previously thought. While local isotropy of the 
r < 10/i _1 Mpc scale correlation function can in princi- 
ple be used as a local geometrical constraint, a standard 
ruler at ^Lgg ~ 130 h^ 1 Mpc has the advantages (i) of 
being little affected by peculiar velocities, and (ii) of oc- 
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curring well into the linear regime where evolution within 
a Hubble time is unlikely. 

Moreover, the detection of the Lt,ss/2 peak (cf fig. 6 of 
Tadros fc Efstathiou 19961 fig. 1 of [Sutherland et al. 1999] ; 
fig. 3 of Mo et al. 1992) implies that both peaks might 
cither be signs of high baryon density (Eisenstein 1998; 
Meiksin et al. 1998; Peebles 1999) or of pre-inflationary 
physics ( Martin fc Brandcnbcrgcr 2000| ), enabling con- 
straints to be put on these. For increased confidence in 
this method, more precise low redshift constraints on large 
scale structure features near ij^gs ~ 130 ft -1 Mpc will be 
highly desirable. Results from the 2 Degree Field Galaxy 
Redshift Survey (e.g. Folkes et al. 1999|), the S loan Digital 
Sky Survey galaxy sample (York et al. 2000), and the 6 



Degre e Field galaxy survey (e.g. Mamon 1998 ) may help 
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